Theoretical description of time-resolved photoemission spectroscopy: application to 

pump-probe experiments 
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The theory for time-resolved photoemission spectroscopy as applied to pump-probe experiments is 
developed and solved for the generic case of a strongly correlated material. The formal development 
incorporates all of the nonequilibrium effects of the pump pulse and the finite time width of the 
probe pulse. While our formal development is completely general, in our numerical illustration 
for the Hubbard model, we assume the pump pulse drives the electrons into a nonequilibrium 
configuration, which rapidly thermalizes to create a hot (quasi-equilibrium) electronic system, and 
we then study the effects of windowing that arise from the finite width of the probe pulse. We 
find sharp features in the spectra are broadened, particularly the quasiparticle peak of strongly 
correlated metals at low temperature. 

PACS numbers: 71.27.+a 71.10.Fd 71.30.+h 79.60.-i 



Pump-probe, femto-second (and recently, atto-second) 
time-resolved photoemission spectroscopy (TR-PES), 
and time-resolved angle resolved PES (ARPES) tech- 
niques can directly examine the excited state nonequilib- 
rium dynamics of electrons in solids llj, including some 
strongly correlated electron systems 0, 0] ■ In these ex- 
periments, an intense pulse of radiation "pumps" the sys- 
tem into a highly excited nonequilibrium state. After 
a variable time delay, the system is subject to a weak 
"probe" pulse of higher energy photons, ejecting pho- 
toelectrons which are detected with energy (and angle) 
resolution. 

Conventional (continuous probe beam) ARPES in lay- 
ered materials can be well approximated Q as a di- 
rect measure of the momentum and frequency dependent 
lesser Green's function [j| of the electrons in the layers 
(i.e., their spectral function multiplied by the Fermi func- 
tion). In more isotropic (three-dimensional) materials, 
the spectral function gets averaged over k z , the compo- 
nent of the momentum perpendicular to the layers. The 
theoretical situation is less clear for pump-probe photoe- 
mission. The interpretation generally used is that the 
pump creates "hot electrons" in quasi-equilibrium at a 
higher effective electronic temperature (T e i ) compared to 
the lattice (phonons), which then cool gradually, so that 
the probe PES essentially measures "equilibrium" lesser 
functions at different values of T e [ for different time de- 
lays between the pump and the probe pulses. For ex- 
ample, recent TR-ARPES experiments Q on the lay- 
ered material lT-TaS2 [believed to be an unusual, charge 
density wave (CDW)-induced, Mott insulator,] were in- 
terpreted in this way, using (equilibrium) spectral func- 
tions from a dynamical mean-field theory (DMFT) [f| 
treatment of the 2-d Hubbard model for a range of tem- 



peratures and fillings, chosen as fitting parameters, for 
the different time delays used. Such an approach, while 
reasonable in many contexts, avoids addressing two im- 
portant questions connected with (1) the nonequilibrium 
dynamical aspects of the experiment, and (2) the effects 
arising from the finite width of the probe pulse (in the 
time domain). 

In this letter, we address these questions — the first, 
by developing a precise formulation for what is actually 
measured in the TR-PES experiments, making future cal- 
culations of the nonequilibrium effects possible; and the 
second, by studying the effect of the probe pulse width, 
but using equilibrium spectral functions at high effective 
electronic temperatures as in the previous studies. 

We begin by developing a theoretical treatment for 
what is measured in TR-PES. We assume that the sys- 
tem, modeled by a quantum many-body Hamiltonian TL, 
is in equilibrium at a temperature T before the pump is 
turned on. It is represented in the distant past (t — ► — oo) 
by an ensemble of the (many body) eigenstates |\l/ ra ) 
of TL, present with the Boltzmann probability p n — 
Z~ x exp[— E n /{kBT)\ where E n are the corresponding 
energy eigenvalues, and Z = ^ B exp[-£„/(fer)] is the 
partition function. Turning on the pump pulse modifies 
TL into a time dependent Hamiltonian TL puinp (t) whose 
precise form depends on the way one models the interac- 
tion of the pump radiation [represented by the vector po- 
tential A pump (r, t) whose t dependence includes its turn- 
ing on and off] with the electrons of the system [i| . Then, 
at time to, just before the probe is turned on, the sys- 
tem is represented by the ensemble of states [^(in)) = 
U (to, — oo)|^„) (with the same Boltzmann probability as 

above), where £/(t 2 ,ti) = %{exp[~i J t * 2 dt'TL pnmp (t') / fr]} 
is the unitary time development operator of the system in 
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the presence of the pump radiation. Here T t is the time 
ordering operation that handles the non-commutativity 
between H pump (i') at different times. {|^(*o))} ac t s as 
the ensemble of 'initial' states for the quantum transi- 
tions generated by the time-dependent probe pulse. By 
explicitly incorporating the time-evolution operator of 
the quantum system with the pump, we include all pos- 
sible noncquilibrium dynamics into the formalism. 

Similarly, when the probe is turned on, the Hamilto- 
nian is modified to 7i pU mp(i) +W pro be(i) to appropriately 
include the vector potential A pro t, e (f, t) arising from the 
radiation field of the probe pulse (it is simplest to think 
of the pump and probe pulses as being active at nonover- 
lapping times, but this is actually not a requirement). If 
the system is in a particular 'initial' state |\l/^(to)), then, 
at any later time t it will evolve to the 'final' state 



(1) 



where U(t,t ) = % cxp[-i dt'{H pump (t') + 
7~Lpy:obe(t' )} / fi] is the full time development opera- 
tor in the presence of the probe pulse. When we have 
nonovcrlapping pump and probe pulses, Ti pU m P (t') = Tt 
for t' > to- After the probe pulse has been turned off, 
the probability that a photoelcctron is in a final state 
with momentum k e = k e k e (in a momentum interval dk e 
and solid angle dfl^ ) is given by 



lim 



(k e ) 2 dk e dflt 
(2^~^ 



P(t); P(t) = ^2p n \{^ m ;k e \<(t)) 



(2) 

Here, as appropriate to photoemission from an experi- 
mental sample with a surface, \*$> m ;k e ) is well approxi- 
mated as a direct product of the many-body eigenstate 
l^rn.) of the initial (and final) time-independent Hamil- 
tonian TL, and a one-electron scattering state which is a 
free electron state of momentum hk e outside the sample. 
Since the eigenstate \^m) in which the system is left is 
not determined in the experiment, and the initial state 
can be any one of the ensemble of initial states with prob- 
ability p n , Eq. ^ includes an unconstrained sum over m, 
and a sum over n weighted by p n . 

We presume that the pump pulse is intense enough 
that it needs to be treated non-pcrturbatively, but that 
the probe pulse is weak enough that Ti pm bc(t') can be 
treated by perturbation theory. Hence, to leading order 



U(t,t )~U(t,t )[l- 



dt'u\t',t )n plohe (t')u(t',t )} 

(3) 

The pump radiation is chosen such that its photons do 
not have sufficient energy to overcome the work function 
W of the sample and eject photoelectrons. Hence the 
photoemission process arises only from the second term 
in Eq. ([3]). The component of H. pm bc(t') which causes 
the absorption of a photon of momentum q (frequency 



ujjj = cq ) and the ejection of an electron of momentum 
k inside the system as a photoelectron of momentum k e 
is of the form: 

n pTobe (t') = J dk z M(q,k,k ez )s(t')a ? ci^ (4) 

Here aq is the annihilation operator for a photon with 
momentum q , and for simplicity, we have dropped the 
spin indices of the electron operators. The probe shape 
function, s(t'), captures the time dependence of the probe 
envelope, including its turning on and off. M(q,k,k ez ), 
the matrix element for the process, depends on the details 
of the modeling of the sample, especially its surface [Icj ]. 
This choice assumes the sudden approximation, where 
the photoelectron rapidly moves out of the sample. To 
leading order in 7i pro be (and using the factorization of 
|*m! k e )), we have 



;* m ;* e |^(t)) 



dk z M(q,k,k ez ) / dt' s{t')e 
J t 



-iujt' 



x (* m |C/ t (i',to)c,C/(t',t )|^(to)) • (5) 



Here hXu = hw^ — (hk e ) 2 / (2m e ) — W is the energy of 
the excitation left in the system after the photoemission 
process. For probe photons of fixed direction and en- 
ergy, and for a given material, specifying uj determines 
k e , hence the probability P(t) in Eq. ([2|) is a function 
only of t, u> and k e . Using the properties of the time de- 
velopment operator, and the completeness of {l^m)} it 
is straightforward to show that 



P(t,co,k e ) 

with 

I(t,uj,k e ;k z ,k' z ] 



dk' z / dk z M(q, k, k e 



1 

- w 

x M*(q,k' ' ,k ez )I(t,u,k e ;k z ,k' z ) (6) 



dt" / df's^Os^'K^"" 
x GS-(t',t"). (7) 



Here k = (ku,k z ), k' = (kn,k' z ) , and G < is the well 
known two-time (noncquilibrium) lesser Green's func- 
tion Q given by 

n 

x c u Utf, -oo)\* n ) (8) 
= iZ-^le-Kl^cl^cnit')] (9) 

where ct (t") and c^(t') are the electron creation and de- 
struction operators in the Heisenberg picture appropriate 

to 7~C p um P (t) • 

ct,(t") = U{-oo,t")ciU{t",-oo); (10) 
cg(i') = U(-oo,t')c^U(t', -oo). (11) 
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Frequency tu/t* Frequency w/t* Frequency w/t* 

FIG. 1: (Color online) Photoemission spectra for the Hubbard model at half filling and U = 3 for panels (a) and (b) (metal) 
and (7 = 4.2 for panel (c) (insulator). The temperatures are T = 0.00539 for panel (a) and T = 0.0189 for panels (b) and (c). 
The continuous beam photoemission curve is in black and is labeled by G k (uj), while the pump-probe curves for different probe 
widths are in other colors and labeled by their widths T. Note that there is very little difference between the continuous beam 
and r = 10 curves in panels (b) and (c). The insets focus on the low-energy structure, where there are differences. 



In the pumped case, this needs to be calculated using 
nonequilibrium Keldysh contour-ordered Green's func- 
tion techniques § on the Kadanoff-Baym-Keldysh con- 
tour in the complex time plane. Of course, if the pump 
pulse and the envelope function s(t) vanish for t > t\, 
then P{t) becomes independent of time for t > t\. 

In the equilibrium case, when s(t) = 1, GS(t',t") is 
only a function of (t! — t"). Furthermore, in a highly 
anisotropic layered system, its k z dependence can be ne- 
glected. Then the ARPES transition rate (transition 
probability per unit time, which is what is relevant as 
the continuous beam probe pulse photocmits electrons 
at all times at a constant rate) is proportional to 



I(t,uj,k e ) 

t— >oot -^-oo [t — to) 



lim lim 



-iG$(w)=A i (u,)f(w) (12) 



which is the standard result. 

In our numerical results, we focus on a material that 
can be well approximated by the d — > oo Hubbard model, 
for which the DMFT is exact IfJ. The Hamiltonian is 



H = - 



t* v t 



(13) 



Here electrons hop between nearest neighbor sites on 
an infinite dimensional hypercubic lattice (leading to a 
Gaussian noninteracting density of states) and they have 
an opposite spin repulsive interaction U when two elec- 
trons are on the same site. In the rest of our paper, all 
energies are in units of t* . We consider only TR-PES, 



where we integrate over the direction of the photoemit- 
ted electron. For simplicity we neglect the dependence 
of the matrix element M on momenta [l(| and treat it 
as a constant, and we ignore all k z dependence in /. We 
also assume that the pump pulse acts solely to heat the 
electronic system which rapidly thermalizes into a quasi- 
cquilibrium distribution at an effective electron temper- 
ature T e i, so that we can focus on the effect of the win- 
dowing due to the probe pulse envelope s(t). 

We solve the DMFT problem using the numerical 
renormalization group Q for the retarded Green's func- 
tion. We take A = 1.8 and keep 800 states per iteration 
on the Wilson chain. Once the retarded Green's function 
is found, we multiply the imaginary part by —2if(w) to 
get the lesser Green's function, then we Fourier transform 
to real time to find G < (< — t'). Next we evaluate Eq. (J7J), 
which is proportional to the TR-PES spectra given the 
above assumptions. We take the shape function s(t) to 
be a Gaussian s(t) = exp[— (t — t) 2 /T 2 ]/(T^/t:) with a 
varying width T and a peak time [ll| i; We work at half 
filling, with various U values to examine both metallic 
and insulating phases. Our PES signals are normalized 
so that the integrated weight in all spectra are identical. 

In Fig.Q] we compare the continuous beam PES [s(i) = 
1] with the pump-probe PES for Gaussian probe shape 
functions of varying width. Panels (a) and (b), with 
U = 3, correspond to a strongly correlated metal, and 
panel (c), with U = 4.2, to a Mott insulator. While the 
broader higher energy features of the spectral function 
are determined accurately once the width is on the or- 
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Frequency w/t* 




Frequency cj/t* 



FIG. 2: (Color online) Time-resolved photoemission spec- 
tra for the Hubbard model at half filling and (7 = 3 (top, 
width is 5) or U = 4.2 (bottom, width is 2) at temperatures 
T — 0.00539 (dashed line) versus the equilibrium continuous 
beam PES at different Ts. Note how the windowing effect is 
different from thermal broadening effects, especially for pos- 
itive energies. The two low-temperature curves in panel (b) 
are almost overlapping except at positive frequencies, where 
there is a small difference. 



der of r = 5, sharp features like the quasiparticle peak, 
or kink-like features at higher temperatures (see inset at 
frequencies near ±0.2), require much broader windows 
to be accurately captured. We see similar results for the 
Mott insulating phase, where the spectra have no sharp 
features, except for some positive frequency features due 
to the thermal excitations at high T as in the inset; re- 
sults at lower T are essentially the same, but have no 
thermally activated peaks at positive frequency. 

It is interesting to ask whether the broadening effect 
of the finite probe pulse width can be mimicked by a fi- 
nite temperature thermal broadening. We examine this 
in Fig. [21 where we compare the metallic PES (top panel) 
for T = 0.00539 and a probe pulse width of T = 5 



(dashed line) with the continuous beam PES for various 
temperatures. We see the spectral shapes are quite dif- 
ferent. The reason is that the high temperatures needed 
to broaden the quasiparticle peak to the same extent as 
obtained from the windowing effect of the probe pulse 
also generate higher energy upper Hubbard band contri- 
butions. The windowing effect alone just does not have 
these. This is further supported by the results for the 
insulating phase (lower panel), where we see similar ef- 
fects. 

In conclusion, we have developed a full many-body the- 
ory under the conventional photoemission assumptions 
for pump-probe time resolved photoemission that takes 
into account two new effects: (1) the noncquilibrium dy- 
namics brought on by the large fields in the pump pulse 
and (2) the windowing effect due to the finite width of 
the probe pulse. We find that the PES signal including 
all the nonequilibrium effects can be represented in terms 
of integrals of the lesser Green's function in the presence 
of the pump pulse, and that the windowing effect can 
make it difficult to extract narrow energy features unless 
the pulse width is wide enough in the time domain. This 
latter effect is quite different from the effect of raising 
the temperature, which also causes broadening, but of a 
qualitatively different character. Given the experimental 
parameters and lowest-energy bandwidth of lT-TaS2, it 
is possible that the windowing effect could be playing a 
role in that system. 
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